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Recent measurements of the tt^tt invariant mass distribution at RHIC show a shifted peak for 
{3J[J' the p meson in IQQA GeV in peripheral An + Au and even in p + p collisions. A recent theoretical 

^ ^ study based on a picture of in-medium production rates of pions, showed that a large shift could 

result from a combination of the Boltzmann factor and the coUisional broadening of the p. Here 
we argue that the two-pion density of states is the appropriate quantity if one assumes a sudden 
0^ \ break-up of the system. Methods for calculating the density of states which include Bose effects are 

04 , derived. The resulting invariant mass distributions are significantly enhanced at lower masses and 

the p peak is shifted downward by ~ 35 MeV. 
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, One of the most compelling motivations for studying heavy ion collisions is the prospect for observing the restoration 
of chiral symmetry. The spontaneous breaking of chiral symmetry is accompanied by the creation of a quark-antiquark 
condensate whose coupling to nucleons is responsible for the great bulk of the nucleon mass, and is therefore responsible 
for most of the mass of the universe. The transient nature of the heavy-ion reaction precludes a detailed investigation 
of all the quasi-particle modes in the highly excited collision volume. However, the p meson is unique for it typically 
^ decays inside the spatial region where the vacuum structure might undergo novel changes. A neutral p decays with 
99% probability into a tt+tt^ pair and decays with a small probability into e"*"e^ or p.'^pi^ pairs. The electromagnetic 
channels are especially useful because dilepton pairs will largely leave the collision volume unscathed by interactions 
with the thousands of other constituents. Since the p has the same quantum numbers as the photon, the invariant 
mass spectrum of dileptons is dominated by the p for masses between 600 and 800 MeV. Experiments at the CERN 
^ SPS for e"'"e~ and pC^ n~ 0,0 suggest that the p has either dissolved "l] (as would be expected in a quark-gluon 
■ plasma), has moved down a few hundred MeV (due to chiral symmetry restoration), or has been broadened via 
collisions by many hundreds of MeV lS)] • 

Recently, the possibility of studying in- medium properties of the p meson through the tt+tt" channel has been 
discussed [3, 111 • Unlike dileptons, pions are not penetrating probes and are likely to re-interact before they escape. 
Since temperatures fall to near 100 MeV at breakup, where the p/ir ratio falls to a few percent, the chance that a 7r+ 
is accompanied by a ■k~ that originated from the same p, rather than a charged pion from a different source, is only 
a few percent. Thus, a background subtracted invariant mass distribution should have a p peak that comprises only 
a few percent of the integrated distribution. 

The STAR collaboration at RHIC has measured such a peak in pp collisions, and for the first time, in peripheral 
relativistic heavy ion collisions 9]. A surprisingly significant downward shift of the mass was observed even in pp 
collisions, especially at low pt, and an even larger shift was observed in peripheral Au -I- Au collisions. Results are not 
yet available for central collisions where it is more difficult to observe the peak since the p/n ratio falls. Eventually, 
the p peak should also be measured for central collisions given sufficient statistics. 

In references and the mass distribution was predicted by considering the in-medium rate of p decays into 
7r+7r~ pairs, dN/dMcPxdt. This is the same approach as has been applied for dilepton studies. In reference 0, these 
rates were corrected for collision broadening and for Bose effects. Collision broadening was shown to be particularly 
important in moving strength to lower-lying masses. However, emission of pions is of a fundamentally different 
character than that of dileptons. First, the final-state distribution is not necessarily proportional to the decay rate 
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since the decay rate is often balanced by a formation rate of similar magnitude. Secondly, collisional broadening can 
not be applied in the same manner since measurements are made in the asymptotic state. Finally, the presence of the 
p alters the two-pion scattering partial waves at non-zero separations which should affect the mass distribution. As we 
will demonstrate, the production-rate calculations of and provide different results than a freeze-out prescription 
which is governed by the available phase space. 

If the last strong interactions felt by the two pions used in the distribution can be considered sufficiently hard 
to statistically sample the outgoing phase space, the two-pion density of states should govern the invariant mass 
distribution. Other non-randomizing interactions, identical particle symmetrization and mean field interactions, 
would then serve to modify the density of states. Although these ps probably decayed during the breakup stage, 
which is well below the critical density, the decaying p mesons might still sample a region where mean-field effects, 
i.e., in-medium mass shifts, are not negligible. Since pions are Goldstone bosons, they probably leave the region 
with their energy and momenta unchanged during their exiting trajectory, and one expects that a modification of the 
2-pion invariant mass distribution would reflect the in-medium modifications of the p rather than those of the pion. 
Whereas, this mass shift may be on the order of 100 MeV at high temperature, it is unlikely to be much more than 
25 MeV at breakup when densities have fallen well below nuclear density. 

It is not the aim of the current paper to model the in-medium mass shift of the p, but rather to investigate how the 
invariant-mass distribution would look in a thermal description based entirely on the two-pion density of states, and 
the associated Bose effects. In the next section, methods for calculating the two-particle density of states are presented 
along with a comparison with the functional forms one would expect from rate calculations. After convoluting with 
the Boltzmann weighting, we find that the p peak is shifted downward by ~ 30 MeV relative to the nominal p mass. 
The shift is due to three factors, the Boltzmann weighting T^, the fact that the density of states peaks below the 
p mass, and the inclusion of other partial waves. Bose-Einstein effects also enhance the distribution at lower masses 
Hi llll ll^ , especially for heavy ion collisions where the pionic phase space filling factors are approaching unity [isL ll^ . 
In section UTTl methods are presented for including Bose effects into the two-pion density of states. The resulting mass 
distribution is strengthened at lower invariant masses, but the peak did not shift appreciably. 



II. INVARIANT MASS DISTRIBUTIONS FROM THE TWO-PION DENSITY OF STATES 

Since the first measurements of the p meson |0, 0|, the masses and widths have fluctuated by several MeV 
depending on the analysis. Currently, the Particle Data Group assigns a nominal mass of 771.1 MeV and a width 
of 149.2 MeV with uncertainties for each number being near 1 MeV. The p mass has been determined from a 
number of means, e^e~ — > tt+tt" reactions, pp collisions, and up — > ttttj* reactions p^ . Electro-production of the p is 
complicated by the interference with the cj ^ 27r channel which constructively interferes with the p^ channel since 
the electromagnetic coupling violates isospin conservation. Since pp collisions are typically highly inelastic, extracting 
the p mass is complicated by the same factors that complicate the study in a heavy-ion environment. In the irp ■n'Kp 
reaction, the proton is treated as a source of pions which are assumed to scatter elastically with the incoming pions. 
In fact, 7r+7r~ phase shift analyses have been successfully performed. The cross section for tt+tt" reactions should 
have a Breit-Wigner form, 



where q is the momentum of either pion in the center-of-mass frame (M = 2 + 9^ ) , and Ylj is the imaginary part 
of the one-loop self energy of the relativistic propagator. 
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n/,0 = ToAfo. (3) 

(4) 

Here, Fq is the nominal width, and go is the momentum required to provide the nominal mass, Mg. The last term, 
P{q/Qo) is a form factor whose exact form is in doubt [TH . The product q^a peaks precisely at the nominal mass 
irrespective of the form factor. 
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FIG. 1: The spectral function of the p (dotted line) is broader than the Breit-Wigner form (dashed line). The total density of 
states, including effects of modifying the outgoing partial waves, is noticeably shifted to the left relative to the other forms. The 
difference is especially noticeable at small invariant masses, where the three forms rise proportional to , (f' and q respectively. 



The spectral function of the p is related to the imaginary part of the propagator. 



2M^ 1 

TT (A'P - M^) + ail 

2M 



S,{M) = (5) 



BW{M) (6) 



Here, Sp is usually associated with the number of states available to the p with a given mass. The real part of 11 is 
being ignored for the current discussion. Since the Breit-Wigner function, BW{M), always peaks at M = Afg, and 
since XIj/M is rapidly growing with M near the p mass, the spectral function always peaks below the p mass. Setting 
the form factor in Eq. (|2Jl to unity, the peak of the p spectral function shifts downward by 5 MeV. Applying some of 
the different expressions for 11/ discussed in 11] may result in the peak being shifted further downward, perhaps as 
much as an additional 5 MeV. 

The change in the total density of states can be expressed in terms of phase shifts, [2(1 i2H : 

Given the relation between the phase shift and the self energy, one can express Ap in terms of the self energy. 

The first term is the spectral function of the p, which is often associated with the probability of having a p meson 
of mass M . Together, the two terms describe the entire correction to the density of states, including the effects of 
modifying the outgoing partial waves. 

Figure^illustrates the importance of using the correct expression for the density of states. The spectral function of 
the p is peaked below the Breit-Wigner function, and the total density of states is peaked even lower. The difference 
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FIG. 2: The thermal mass distribution of ■k'^-k~ pairs is shown for three calculations, both at T = 110 MeV and T = 170 
MeV. Using the full density of states as calculated from taking derivatives of phase shifts results in a broader distribution for 
the i = \ channel (full line) than using the p spectral function (dotted line). Including all s, p and d channels (dashed line) 
provides significant strength at low invariant masses due to the s wave channels and moderate strength at higher masses from 
the d wave channels. 



is especially strong at low invariant masses, as the Breit-Wigner function rises as g^, the p spectral function rises as 
(f" and the pionic density of states rises as q. This relative scaling with q would hold for any p-wave interaction. 

Thus far, the distribution of masses has not incorporated the Boltzmann factor, which should push the peak even 
lower with the thermal weight, e~*^/"^ |10|. More precisely, one needs to integrate over the modes in momentum space 
due to relativistic effects, 

This should represent the background-subtracted 2-pion invariant-mass distribution. As can be seen in Fig. |21 the 
Boltzmann weight pushes the distribution increasingly downward for lower temperatures. The upper panel shows the 
mass distribution assuming a temperature of 170 MeV, which is a reasonable temperature for thermal models of 
collisions, while the lower panel shows the result for a temperature of 110 MeV, which may be reasonable for the 
breakup temperature in central heavy ion collisions. Calculations using both the p spectral function and the two-pion 
density of states are displayed to illustrate the importance of choosing the appropriate form for the density of states. 
The Boltzmann factor greatly magnifies the enhancements at low Af , to the point that a second peak appears for 
lower temperature. 

The TT+TT^ density of states is also affected by phase shifts in other channels. For tt+tt^, the s-wave channel is 
split into two isospin components, 2/3 weight for / = and 1/3 weight for 1 — 2. The / = channel is particularly 
important as it corresponds to the mythical a meson. Although phase shift analyses do not reveal a sharp peak as in a 
resonance l23ll23. 125ll26l| . the phase shifts are considerable, rising steadily from zero at threshold to approximately 
90 degrees at M = IMk ~ 1 GeV, where the kaon channel opens. At the two-kaon threshold, the behavior of the 
phase shifts becomes complicated and an inelastic treatment becomes warranted. Since one uses derivatives of the 
phase shifts to find the density of states, interpolating data for phase shifts can be dangerous due to noise in the 
experimentally determined phase shifts. Thus, we apply a simple form that describes the general behavior. 



(5/=o,s=o = ag + h{M - 2m^). 



(12) 
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The first coefficient a is tlie scattering length, which is small due to constraints from chiral symmetry. The number 
varies throughout the literature by several tens of percent. We use the value, a = 0.204/m^ [23|. The second term 
does not contribute to the scattering length, as (M — 2to^) ~_2^ at low q. Choosing h = 9.1 x 10"'* GeV~^ crudely 
reproduces experimental phase shifts, which are reviewed in |22j . Since these phase shifts rise half as far as those in 
the delta channel, have one third the spin degeneracy, and have a 2/3 weight in the 7r"'"7r~ channel, they are noticeably 
less important than the p channel in affecting the overall density of states, unless one is near the two-pion threshold 
where p-wave interactions vanish. 

Other phase shifts also contribute: (/ = 2,^ = 0), (/ = 0,^ = 2) and [I = 2,1 ~ 2). Since none of these phase 
shifts exceed more than a few degrees, they make nearly ne glig ible contributions to the density of states. For the 
(/ = 2, = 0) channel, we apply an effective range expansion |28| . 

cot 5 = — + ]-Rq, (13) 
qa 2 

where a — —0.13 MeV^^ and R = 1.0 MeV~*. The d wave is also composed of / = and / = 2 pieces. For the 
(/ = 0, ^ = 2) piece, the data are rough, and we make a simple expansion, 

(5/=o,fc2=cg^ (14) 

where c = 6.2 GeV~^. The parameter a is uncertain to the 50% level. For the {I = 2J, — 2) partial wave, we use an 
expansion ^2^, 

5i=24=2 = -8.4g'5 + 12. GeV"^ (15) 

None of the these three channels are well understood, but none have a substantial impact at or below the p region of 
invariant mass. 

FigureElalso shows the invariant mass distribution of a thermal ensemble with T — 110 MeV using all the s, p and 
d channels. The s-wave contributions are non-negligible near the p mass, and dominate near the two-pion threshold. 
The d-wave contributions matter only for masses near or greater than 1.0 GeV. 



III. BOSE EINSTEIN CORRECTIONS 



Bose Einstein corrections should preferentially enhance low-mass pairs since low-mass pairs are more likely to 
include a low-momentum pion. This has been investigated within the context of the p peak as well as the influence 
on Z boson decay modes |lll |. In this section, we present a means to include Bose enhancement effects which are 
consistent with the statistical picture described in the previous section. 

In order to demonstrate Bose enhancement effects, we revert to the fundamental definition of the two-particle 
density of states. 

J [V \ )" (17, 

27r ^AI-Ho + ie\M-H„ + ieJ ^ ' 

We will work in the two-pion rest frame, so the trace would cover all two-pion states that have total momentum 
zero. When including Bose effects, one would sum all such two-pion states, plus average over the distribution of other 
identical particles whose probability of being populated is 

/(q) - T^^rh. - + (1^) 
-L - /oiqj 

Thus, (1 -I- /) can be considered as an Bose enhancement factor while /o is the phase space filling factor if Bose 
statistics were neglected. 

Using the cyclic property of the trace, Eq. (|16|l can be written in terms of a derivative with respect to M, 

p(M) = «o(M) + — 3— Tr V i { ^ ^ (19) 

' ' 2tt dM ^n\M-Ho + iel ^ ' 



n—1 
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One can expand the n terms and note that the sum includes aU possible orderings of n factors where each factor is 
either the principal value piece, which is real, or the imaginary part which is proportional to the density of states. 
One could restrict this sum to cover only those terms where the first factor is the imaginary part and multiply by a 
factor of n/Np, where Np is the number of times that i-KS{M — Hq) appears in the term. The sum over n can then be 
transformed into a sum of all possible numbers of appearances of the real part. 

p{M) = po(M) + J^TV Wp ^'""^^^ - 
Here, TZ is often referred to as the i?-matrix. This can be written in terms of a logarithm, 

. 1^ , / 1 + ITTdiM - Ho)'Jl\ , , 

Thus, the density of states is determined completely by a single matrix, 

r = 7rpoiM)Tl, (24) 

which is evaluated only for those states whose energy equals M. In a partial-wave basis, r is related to the phase 
shift, T — tan (5. In a plane-wave basis, the matrix r links one direction of the relative momentum with another, i.e., 
the matrix should be written with indices, tq^^q^. 

The presence of other particles alters r. Each matrix element V used to construct r is modified by the presence of 
other particles by the Bose enhancement factor. 



nqi,-qi;q2,-q2) ^ ^qi, -qi; q2, -q2) V(i + /(qi))(i + /(-qi))(i + /(q2))(i + /(-q2)). (25) 

If the intermediate states contained in the definition of TZ are not affected by the phase space density, one can scale 
T in the same manner as V. Then, given the fact that each state appears in both the bra and ket, one can modify r 
in a simple manner to account for Bose effects, 

T(qi, -qi; q2, -q2) = To(qi, -qi; q2, -q2)(l + /(qi))(l + /(-qi))- (26) 
The density of states is comprised of integrals of a cyclic nature, 7„, 

n=l,3,5--- 

/„(M) = /^^...^ro(Oi,02)ro(fi2,fi3)---ro(n„,Oi) (28) 
(1 + /(qi))(l + /(-qi)) •••(! + /(q„))(l + /(-q„)). 

Unless the momentum of the pair P = 0, the phase space densities will be sensitive to the direction of the relative 

momentum il. 

For a purely s-wave interaction, tq has no angular dependence and /„ easily incorporates Bose effects, 

In=(ro J ^(l + /(q))(l + /(-q)))". (29) 
The correction to the density of states is then 

A.(M) = I'i^. (30) 

r = tan^l ^(1 + /(q))(l + /(-q)), (31) 

where 5 is the phase shift as measured in the absence of Bose modifications. 
For a p wave interaction, tq has the angular dependence, 

To (Hi, 1^2) = foQi ■ h- (32) 
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By choosing a coordinate system where the z axis is parallel to the total pair momentum, there is reflection symmetry 
about the x, y and z planes. By making use of the identity, 

— {A-h)(b-C)F{nb) ^ A'-C (33) 
Ait 

K = AF,, (34) 



J ^F(17)cos2 0sin2 0, (35) 
—F{n)sm'^(t)sm^e, (36) 
= 1 ^Fin)cos'e, (37) 



one can iteratively perform the integral in Eq. (|28|l . 



In{M) = (foi^.)" + (foi^,)" + (foi^.)" , (38) 

(39) 

Using F{Q) = (1 + /(q))(l + /(-q)), one can calculate Ap, 

Ap(M) . i f ^ + ^ + ^) , (40) 

— (1 + /(q))(l + /(-q))3sin2 c^-sin^ 9, (41) 

— (1 + /(q))(l + /(-q))3cos2 0sin2 9, (42) 

= tantJ y ^(1 + /(q))(l + /(-q))3cos2 0. (43) 

The calculation of Ap{M) must be repeated for each value of the total momentum since /(q), which is defined in the 
two-pion rest frame, changes when the total momentum is changed. 

The p-wave and s-wave corrections to the density of states do not interfere with one another since they have 
opposite parities and (1 + /(q))(l + /(— q)) has even parity. However, calculation of the £ — 2 contributions would 
be complicated by the fact that the elliptical distortion of the Bose enhancement factors would mix the £ — 
and £ = 2 contributions. For the calculations here, the £ = 2 contributions were calculated by assuming that the 
Bose enhancement factors were independent of il, then using enhancement factors which had been averaged over all 
directions of fl. 

The mean Bose enhancement, ((1 + /i)(l + /2)), is shown as a function of the invariant mass and momentum of the 
decaying in Fig. O assuming a breakup temperature of 110 MeV and an effective chemical potential of 90 MeV. 
The enhancement has been averaged over the directions of the relative momentum. The enhancement is largest for 
low-momentum, low-mass pairs since these pions most strongly sample the region of high phase space density. For 
higher invariant masses, the Bose enhancement is actually stronger for higher pair momenta, as it allows one of the 
outgoing pions to have low pt and sample the high phase-space density region. From viewing Fig. O it is clear that 
the Bose modifications to the invariant mass distribution would be more acute if experiments were to focus on pion 
pairs with low total momentum. 

Bose corrected densities of states are shown in Fig. 0|for T = 110 MeV and p — 90 MeV. A non-zero chemical 
potential was used to account for the relative overpopulation of pionic phase space which may result from rapid 
cooling ^33| and might be magnified by the effects of chiral symmetry restoration [sj ■ Analyses of nn correlations 
from RHIC indeed point to high phase space densities [isl ll^. especially for central collisions of heavy ions. As 
expected, lower-mass states were more enhanced by Bose effects. Since the density of states was proportional to the 
derivative of (tan5((l -I- /)(1 + /')), and since the averaged phase space filling factors generally fall as M increases, 
the density of states was less enhanced for intermediate masses as compared to the no-Bose case. The peak of the 
distribution shifted downward by only one MeV after the inclusion of Bose effects. 

Although the position of the peak was not much affected by Bose effects shown in Fig. 0] Bose effects led to a 
near doubling of the distribution at low masses. These effects are most important at low pt where the phase space 
densities are higher. In pp collisions, a movement of the p peak was observed for low pt pairs|3,|32 which is suggestive 
of Bose effects. However, at the p peak each pion has a relative momentum of ~ 300 MeV/c and will largely sample 
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FIG. 3: The mean values of {(1 + /i)(l + /2)) are shown as a function of the invariant mass and total momentum of the 
outgoing pion pair. The values have been averaged over all directions of the relative momentum. The enhancement factors 
exceed 2.0 for a low values of the pair momentum and invariant mass. 
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FIG. 4: The mass distribution is shown with (solid line) and without (dashed line) Bose effects for T — 110 MeV, ^ = 90 MeV. 
Bose effects enhance the probability of producing low invariant-mass pairs since they are more likely to have low momenta and 
stronger Bose enhancement factors. 



phase-space regions with moderate to low phase space densities. Although the invariant mass distribution is mainly 
altered at invariant masses below the p peak, Bose effects should contribute to washing out the peak by increasing 
the declining background to the peak structure in Fig. ^ 
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IV. SUMMARY 

Our principal finding is that the p peak in the tt+tt^ invariant mass distribution should be approximately 35 MeV 
lower than the nominal p mass, if one accepts the scenario of a sudden breakup that thermally samples the two- 
pion density of states. The shift was the result of convoluting the density of states which is shifted by ^ 10 MeV 
below the p mass with the Boltzmann factor. Given the extra cooling inherent to heavy ion collisions, the breakup 
temperature is probably near 110 MeV, well below the characteristic temperatures used to describe pp collisions. This 
low temperature is responsible for the additional downward shift of the peak in heavy ion collisions. In addition to 
the shift of the peak, the distribution showed significant additional strength at invariant masses near the two-pion 
threshold. This additional strength hinged on using the correct expressions for the density of states, especially in the 
£ = channels. Although the position of the peak was not much affected by Bose effects, Bose effects led to a near 
doubling of the distribution at low masses. 

The thermal model presented here rests critically on a pair of assumptions. First, we have assumed that the breakup 
is sudden, i.e., the last strong interaction experienced by the particles samples the outgoing two-particle phase space. 
Indeed, interferometric measurements do suggest a sudden breakup [iJ, |33, yl] ■ If emission were gradual, e.g., surface 
evaporation, this picture would be invalid. An appropriate treatment of the surface would include the dynamics of 
surface penetration and absorption and might include coUisional broadening. For instance, spectral lines in stars 
are affected by coUisional broadening. The "truth" of the breakup at RHIC probably has elements of both volume- 
like breakup and surface-like evaporation. Thus, the effect of collisions, which played a pivotal role in moving the 
distribution downward in @, requires more study. 

The second assumption inherent to these calculations is related to the neglect of finite-size effects. The enhancement 
factors apphed to small-angle correlation studies are usually based on the outgoing wave function, |(/)(q, r)p [sslls^ . 
For large sources, there is a straight-forward correspondence between the integrated wave functions and the phase 

shifts ila, m, 

i5:(2^+l)g = |^/dV(|<A(q,r)p-|<^o(q,r)p). (44) 

i 

Since the modification of the wave function is mostly confined to a region where qR < tt, one expects that the 
treatments shown here should work well for q > 100 MeV/c, or for masses greater than 400 MeV. For masses near 
threshold, a different approach, based on the actual scattered wave functions, would be warranted. Such an approach 
could also incorporate the effects of the Coulomb interaction between pions. 

Finally, it should be emphasized that other correlations, besides those resulting from the change in the two-pion 
density of states, will play a role in any experimental measurement. Experimental analyses are typically based on 
a like-sign subtraction. This should eliminate global correlations such as elliptic flow which correlate same-sign and 
opposite-sign pairs equally. However, any correlation based on charge conservation should survive the subtraction 
|33.l40| . For every 7r+, there is a ~ 75% chance that local charge conservation will result in an extra tt" being emitted 
with a similar rapidity. This should provide a bump in the like-sign-subtracted invariant-mass distribution that peaks 
for masses near 400 MeV. The ratio of the p peak in the like-sign subtracted distribution to the bump from charge 
correlation is approximately determined by the chance that a given 7r+ had its last interaction with other hadrons 
through the decay of a p". This ratio should be smaller for central collisions since the breakup temperature is lower 
which reduces the p/ir ratio. 

The in-medium mass of the p might be altered by ~20 MeV at breakup. Given that this peak is also spread out and 
distorted as shown in the calculations presented here, it is certainly challenging to isolate the contribution from the 
p and to quote a peak height to a better accuracy than 20 MeV. Upcoming runs at RHIC may increase the statistics 
by more than an order of magnitude. Thus, we believe that there remains a good chance that the p can be studied 
in detail, even in the central collision Au-|-Au environment. 

Finally, we compare the experimentally observed mass shifts to results of our model. In reference 'Oj it was reported 
that the p shift downward in pp collisions by ~ 20 MeV at higher pt while shifting downward by ~ AbMeV at low 
Pt- The shift appeared to be 5 to 10 MeV larger for high-muliplicity pp collisions and perhaps another 5 MeV lower 
for peripheral Au-I-Au reactions. Similar behavior for pp collisions had been reported for y/s = 27 GeV pp collisions 
[s^. The shifts that we extracted were as large 35 MeV, but these calculations assumed a lower temperature, 110 
MeV, and a higher effective chemical potential, 90 MeV, than would be appropriate for pp phenomenology. For a 
temperature of 170 MeV, and zero chemical potential, the shift was in the range of 20 MeV, a somewhat smaller shift 
than what was observed by STAR. It appears that the experimental mass shift is 10 to 20 MeV stronger than what we 
would expect from our approach. But, before this discrepancy can be attributed to novel in-media phenomena, i.e., 
a mass shift of the p, it should be stressed that systematic uncertainties described in |3| are of the order of 10 MeV. 
This problem would be served well by both a higher statistics experimental analysis and a more detailed theoretical 
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modeling. An improved calculation would consider finite-size effects, the influence of other resonances, and the effects 
of experimental acceptances and efficiencies. 
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